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Abstract 

We suggest a new cosmological scenario which naturally guarantees the smallness 
C<") ' of scalar masses and VEVs, without invoking supersymmetry or any other (non- 

o 



gravitationaly coupled) new physics at low energies. In our framework, the scalar 
masses undergo discrete jumps due to nucleation of closed branes during (eternal) 
inflation. The crucial point is that the step size of variation decreases in the direction 
of decreasing scalar mass. This scenario yields exponentially large domains with a 
distribution of scalar masses, which is sharply peaked around a hierarchically small 
value of the mass. This value is the " attractor point" of the cosmological evolution. 



1 General idea 



X 

The radiative instability of scalar masses is the key point of the gauge hierarchy 
problem. In the effective 4D field theory, the scalar masses are quadratically sensitive 
to the ultraviolet cutoff. The only known exceptions to this rule are Goldstone 
bosons. This fact is hard to reconcile with the observed smallness of the weak scale, 
relative to the Planck mass M p ~ 10 19 GeV. So far supersymmetry is the only 
known symmetry that renders masses of elementary scalars radiatively stable. The 
scalar masses are controlled by supersymmetry breaking scale. Given the fact that 
we do not understand the origin of this scale, supersymmetry per se does not really 
explain the origin of the weak scale but rather makes the gauge hierarchy technically 
natural. 

In view of the above, it is crucial to explore other possible mechanisms of scalar 
mass stabilization. In the present paper we suggest an alternative mechanism that 
can guarantee zero or very small scalar masses (and VEVs) without invoking super- 
symmetry or any other non-gravitationaly-coupled new physics at low energies. 

In our scenario, a small scalar mass is selected with probability one during the 
cosmological evolution. This selection works as follows. We construct a simple 



framework in which scalar masses (and VEVs) undergo discrete variations due to 
nucleation of closed domain wall bubbles (branes) during inflation. Values of the 
scalar mass on different sides of the wall differ by a finite step. The bubbles expand 
exponentially fast and create domains of a new vacuum with a new value of the scalar 
mass. New bubbles are created within the old, and the scalar mass changes further. 
Since inflation is known to be eternal PUI2|j the process of wall nucleation continues 
forever, populating the Universe with exponentially large domains having different 
values of the scalar mass. However, not all the values of the scalar mass (VEV) are 
equally probable. In our model, in the absence of gravitational back-reaction, the 
probability is sharply peaked around zero, because the step A0 decreases towards 
small values of the VEV faster than itself. That is, 



where n > is some power. As a result, the density of states diverges for small 
VEV (mass) of 0. 

Thus, in the first approximation, the probability distribution for has an in- 
finitely sharp peak at = 0. We will show, however, that infrared effects, such as 
the Gibbons-Hawking temperature and quantum fluctuation of during inflation, 
can shift the most probable value of the scalar mass (and VEV) away from zero to 
a small value and round off the maximum of the peak. 

2 Cosmic attractors 

To introduce our mechanism, we use a simple toy model. The main ingredients 
are: (1) a scalar field 0; (2) domain walls (branes) charged under an antisymmetric 
three-form field A^, with the field-strength F a p^ = Fe a /3^s- These objects are 
engaged in the following interrelation. The branes are sources for the three-form 
field. The value of the brane charge is determined by the VEV of 0. The VEV of 
is in turn determined by the three-form field strength F. 

These couplings result in the following dynamics. Nucleation of a closed brane 
changes the value of F. The step of change (the brane charge) is determined by 
0. We construct the model so that an increase in F decreases 0, which in turn 
decreases the charge of new branes that can be nucleated. Decrease of the brane 
charge diminishes the minimal step of change in F. As a result, the subsequent 
decrease of requires more steps, and their number diverges towards small values 



Let us discuss this dynamics in more detail. The action of a free three-form field 
in 4D can be written as 



A0/0 oc 0", 



(1) 



Of 0. 




(2) 



It is invariant under gauge transformations 



(3) 
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where B is a two-form. Due to this gauge freedom, F contains no propagating 
degrees of freedom. The solution to the equations of motion is an arbitrary constant 
value of the field strength, 

F = constant. (4) 

The situation changes in the presence of 2-branes, or domain walls, which may 
act as sources for A, due to the following coupling 

q I A, (5) 

where the integral is taken over the 2 + 1-dimensional world-volume and q is the 
brane charge. The role of such branes can be played by the field-theoretic solitonic 
domain walls jSj (see Appendix A), or by fundamental branes of some sort. Their 
precise origin is unimportant for the present discussion. The change of F across the 
wall is given by 

AF = q (6) 

Thus, F can undergo discrete variations due to nucleation of closed branes 

This mechanism can be used to induce spatial variations of the field 0, by cou- 
pling it to the F-form. We do not require that the Lagrangian contains any small 
scale (such as the supersymmetry breaking scale). We allow the (renormalized) po- 
tential of to be the most general function, including all possible interactions with 
the 4-form field F, 

"(*)=(-"•'+ 5 +-) 5 +-)i*i 4 +- p> 

The dimensionless coefficients are not shown explicitly and are assumed to be of 
order one. The couplings linear in F are suppressed by parity symmetry. The 
couplings in (J7J) effectively convert the mass and the VEV of into functions of the 
4-form field strength, e.g., 

2 ~ (m 2 - F 2 /M 2 ) (8) 

We shall assume that the F-independent part of the mass m 2 takes its natural 
value, m 2 ~ Mp. For definiteness, we shall assume the sign of this contribution 
to be negative and the sign of the F-dependent contribution to be positive. Then, 
the F-dependent contribution will lead to a partiall cancellation of the effective 
mass. This mass will take different values in different parts of the Universe, due to 
nucleation of branes charged under F. 

To ensure that the brane charge q is suppressed at small values of 0, we require 
that the system is invariant under a Z 2 n symmetry, which acts on as 

-> e*w 0, (9) 
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and at the same time changes branes into anti-branes and vice versa (leaving the 
3-form A invariant). 1 The coupling of the 3-form to the branes (j5J) should then be 
replaced by 

2+1 A^ A (10) 

Note that for a non-constant the above coupling is not gauge-invariant, and extra 
non-local, terms have to be added to the action to restore the gauge invariance. A 
detailed discussion is given in Appendix A. These additional terms are particularly 
important for understanding of the radiation of A-waves from the branes that occurs 
if the background value of changes in time (see Appendix B). We have shown in 
3J that a non-local action of this type can arise from a local field theory in a two- 
dimensional toy model after integration over some massless fermions. It is not clear 
whether or not this mechanism can be extended to 4D, and more generally, whether 
or not the required type of action can be obtained in effective field theory. We leave 
this question for future investigation. 

In the context of brane nucleation, however, the additional terms in (jlUj) are 
unimportant. In the regime of interest to us here (small 0), each act of brane- 
nucleation changes by a very small amount, so can be regarded as nearly con- 
stant, and additional terms are negligible. 

The magnitude of the F-step between neighboring domains is set by 0, 

A F oc <j> N , (11) 

and so is the change of the VEV of 0, 

A0 2 oc <P N . (12) 

This is the key point of our mechanism. With every step that decreases the VEV 
of 0, we create a region in the Universe where the brane charges are smaller. This 
allows for finer and finer adjustment of the 0-VEV, accompanied by further decrease 
of the brane charges. Thus, the step size of the field F (and therefore of 0) decreases 
and the "level density" grows towards smaller values of 0, and if the process is not 
for some reason terminated, the total number of levels diverges. 

For instance, imagine that we start in a domain where F ~ 1, and 2 = £ < 1 
in Planck units. In these units, the value of the brane charge in that domain is 

qo = i N/2 (13) 

For the sake of definiteness, let us assume that £ ~ 0.1 or so. After the first step of 
brane nucleation, the change in F is 

AF = g ~ e 12 (14) 



1 Alternatively, we could require that A — > —A, while branes are unaffected. In fact the system 
is invariant under both assignments, independently. 
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and the VEV of is partially cancelled to 

2 - - e /2 ) (is) 

Then it will take approximately n = l/^W 2 ) -1 steps to cancel 2 to 2 ~ £ 2 . At 
this point the brane charge becomes 

qnew ~ £ N , (16) 

and now it will take n ~ 1/^ N ~ 2 steps to calcel to 

2 ~ e 3 , (17) 

and so on. In general, the number of steps required to cancel 2 to an accuracy £ fc 
is 

(number of steps) ~ £ Nk/2 . (18) 

All the allowed values of near = have nearly identical vacuum energies, 
and the corresponding regions will therefore occupy equal fractions of the volume 
in the post-inflationary universe. The corresponding prior probability for is then 
simply proportional to the density of states, 

V*((j))d(f)ocd(f)/(f) N . (19) 

Thus, regions with zero mass and VEV of are maximally probable. This 
hierarchy attractor provides a dynamical mechanism for explaining a zero mass of 
an interacting scalar without need for supersymmetry. 

We note that although the wall charge vanishes as — > 0, the wall tension 
remains large, a ~ Mp, and in the limit the walls become simply domain walls 
separating degenerate vacua. Nucleation of such walls is suppressed by a huge 
factor [5] ~ exp(—irM 2 /H 2 ), where H is the expansion rate during inflation. This, 
however, does not change our conclusions, since eternal inflation provides unlimited 
time for the distribution (fT9"j) to establish. 

In order to use this " attractor" mechanism for solving the gauge hierarchy prob- 
lem, we have to overcome the fact that the attractor point is at exactly zero mass 
and VEV of 0. In the following section we will show that curvature corrections 
to the potential V(<p) generally shift the attractor point away from zero to a small 
value of 0. 



3 Small Higgs mass from quantum fluctuations 

In the above analysis we have ignored the effects of the gravitational back reaction 
on the Higgs mass. One possible source of this back reaction is a non-minimal 
coupling to the curvature, 

\<P\ 2 R. (20) 
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This will create an additional contribution to the Higgs mass during inflation, 

Am lurvature ~ & ■ (21) 

Even in the absence of such coupling, will get a thermal-type contribution to its 
mass due to de Sitter quantum fluctuations. This effect is analogous to that of a 
thermal bath at temperature Tqh ~ H (Gibbons-Hawking temperature jSj). In the 
domains where the Higgs VEV drops below H, the corresponding contribution to 
the Higgs mass is 

Am 2 GH ~ H 2 . (22) 

(For <p ^> H, the fields interacting with get masses greater than Tqh and do not 
contribute to the Higgs potential.) 

To analyze the effect of these contributions on our attractor mechanism, we shall 
first consider a simplified picture, where the expansion rate H remains nearly con- 
stant during inflation. (This situation is realized in some models of hybrid inflation 
|7J.) Due to the mass corrections (JUJ), which we shall assume to be positive, 
the Higgs VEV during inflation will not be given by (|5|). but will rather be shifted 
to 

inflationary ~ K " IK ~ Am ^ ^ 

where Am 2 = Am 2 urvature + Am 2 GH ~ H 2 . Now, the discussion in the preceeding 
section indicates that will be driven not to the point where its post-inflationary 
VEV (JHJ) vanishes, but rather to the point where its inflationary VEV vanishes. 
That is, most of the space in the Universe will be occupied by domains where 

* inflationary ~ ^ ' 

Now, it should be noted that light scalar fields with masses m ^ H are subject to 
large quantum fluctuations during inflation. Assuming the minimum of the potential 
is at = 0o, the characteristic amplituce of the fluctuations 50 is generally given 
by [V(0o + 50) — V(0o)] ~ H A jH]. In our case, if the mass of is driven to zero, 
then V(4>) is reduced to the quartic term in flZJ), and we have 5(f) ~ H. This means 
that the field can be driven to zero only with an accuracy ~ O(H), 

inflationary ~ K " F 2 /M 2 - 0(H 2 )) ~ H 2 . (24) 

After the end of inflation, the gravitational (~ H 2 ) contribution to the mass van- 
ishes, and the VEV is shifted to 

itoaay ~ - F 2 /M 2 ) ~ H 2 . (25) 

Thus, the post-inflationary Higgs mass will not be exactly zero, but will rather 
be comparable to the inflationary Hubble parameter. Because of the quantum fluc- 
tuations, the infinite peak in the probability distribution (jTHj) will be smeared, and 
the distribution will be nearly flat, 7 ? *(0) ~ const, for |0| ^ H. At larger values of 
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0, the probability suppression is at least as strong as in Eq. (JTHJ). It can be even 
stronger, due to the effect of differential expansion. Larger expectation values of 
correspond to smaller values of F, resulting in a smaller vacuum energy (both 
because the potential (JTj) is more negative, and because the F-field energy density, 
pp ~ F 2 , is smaller). As a result, the inflationary expansion rate is lower, which 
can lead to an exponential suppression at large 

In order to solve the hierarchy problem, we require that the peak of the proba- 
bility distribution is at ~ 1 TeV. Then the Hubble expansion rate during inflation 
must be H ~ 1 TeV, and the corresponding vacuum energy density ~ (10 11 GeV) 4 . 
This is a constraint that our mechanism imposes on the inflationary scenario. 

We now discuss the dynamics of the model in some more detail. Quantum 
fluctuations of occur on length and time scales 51 ~ 5t ~ H~ x . Fluctuations at 
the locations of domain walls cause fluctuations of the wall charge q, which in turn 
cause variation of the F-form in the adjacent domains. This variation propagates in 
the form of waves (see Appendix B), from the walls to the interior of the domains. 
The wavelengths of these waves are stretched by the exponential expansion of the 
universe, and as a result, the form field will vary on an exponentially large scale 
in the domain interiors. Shorter waves, emitted near the end of inflation, have not 
travelled far away from the walls. The sizes of the domains are huge compared to 
the present horizon, and it will take an exponentially long time for the shorter waves 
to propagate well into domain interiors. 

After the end of inflation, the Higgs mass takes its zero-temperature value, and 
the Higgs rolls away to its new minimum. At this point, the charge of the walls and 
the values of F-form in the adjacent regions change, but once again, this change 
propagates in the form of waves and is confined to the neighborhood of the walls. 
Moreover, the hierarchy is not destabilized even in regions affected by the change. 
The change of F triggered by the wall is proportional to the final charge of the wall, 
which in attractor domains is set by today's value of the Higgs field ~ 100 GeV. 
Thus, even after the new value of the field-strength is established, the change with 
respect to the inflationary value will be AF ~ 10~ 17N M 2 . Already for N = 2, the 
corresponding change in is A0 ~ 100 GeV, which does not upset the hierarchy. 
For N > 2, the back reaction on is negligible. 

Let us finally discuss how the above scenario is modified in more generic models 
of inflation, in which the value of H fluctuates during eternal inflation and then 
gradually decreases during a prolonged slow-roll period. To be specific, we shall 
assume that inflation is of the "new" type, with expansion rate H max -C M p at the 
maximum of the potential and H m i n at the end of inflation. Another important 
parameter is the borderline value between the regimes of eternal inflation and 
slow roll. In the course of eternal inflation, the bubble walls will be exposed to 
H in the range H* <^ H ^ H max . Since the bubble nucleation rate is so low, a 
typical geodesic will go through the whole range many times between successive 
nucleations. This suggests that the accuracy with which the attractor mechanism 
can drive (^inflationary to zero cannot exceed ~ H max - The slow roll period is relatively 
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short and therefore affects only the immediate vicinity of the bubble walls. Thus, 
in order to solve the hierarchy problem, we have to require H max ~ 1 TeV. 

4 Conclusions 

It is usually assumed that the solution to the hierarchy problem requires introduction 
of some new physics at low energies. In the present paper we have provided a 
counterexample to this statement. We have suggested a novel cosmological selection 
mechanism in which small scalar masses and VEVs become attractors during the 
cosmological evolution. 

The key idea is that (1) the scalar mass is dynamically promoted to a stochastic 
variable that undergoes discrete jumps during eternal inflation; (2) the size of the 
minimal step is a continuous function of the "jumping" scalar VEV <fi. That is, 
the order parameter in question controls its own steps. As a result the probability 
distribution is sharply peaked around a small value, for which the step vanishes. 
We call such a value an "attractor". The post-inflationary value of the scalar mass 
is determined by the mismatch of masses during and after inflation, due to the 
Gibbons-Hawking temperature Tqh ~ H, and by the magnitude of quantum fluc- 
tuations of 0, also ~ H. Thus, the observed value of the Higgs mass in our scenario 
is determined not by ultraviolet physics, but rather by the maximal inflationary ex- 
pansion rate H max . Our scenario requires that H max ~ 1 TeV. It is obvious that this 
constraint is not in any respect analogous to the conventional approaches invoking 
supersymmetry or some other new physics around TeV. 

Our model involves non-local couplings of the form field F, and we have empha- 
sized that it is not presently clear how this kind of coupling can be obtained in the 
low-energy effective theory. We note however that this coupling does not seem to 
have any of the pathologies usually associated with non-local interactions. Another 
unusual feature of our model is that it has solutions describing waves of the field 
F propagating at the speed of light. This suggests the presence of massless degrees 
of freedom, which can potentially lead to testable predictions. These issues need 
further investigation. 

Finally, let us note that the "attractor" mechanism can be combined with the 
usual anthropic approach to solving the cosmological constant problem ^Uj, if one 
is willing to introduced more three-form fields. The main technical problem in this 
approach has been to guarantee a sufficiently small step of vacuum energy variation. 
With our attractor mechanism this is trivially achieved, as it generates branes with 
tiny charges. Anthropic selection can also be used to explain the observed Higgs 
VEV in models of inflation with H max 3> 1 TeV. 

Acknowledgments 

We would like to thank Gregory Gabadadze, Jaume Garriga, and Andrei Gruzi- 



8 



nov for useful discussions. The work of G.D. is supported in part by a David and 
Lucile Packard Foundation Fellowship for Science and Engineering, by Alfred P. 
Sloan foundation fellowship and by NSF grant PHY-0070787. The work of A.V. is 
supported in part by the NSF. 

A Branes with charges suppressed by symmetries 

In this Appendix we shall give a more detailed discussion of branes and domain 
walls with field-dependent three-form charges. Below Planck energies, such branes 
can be treated as fundamental objects, or be explicitly constructed as field theoretic 
solitonic domain walls, as it was suggested in [3J. The dynamics of three- form 
fields in 3 + 1 dimensions is in many respects analogous to the (1 + l)-dimensional 
electrodynamics, with electrically charged particles playing the role of "branes". So 
we shall first review our mechanism in a simplified 2-dimensional example, and then 
generalize to (3 + 1) dimensions. For simplicity, we shall work in Planck units, and 
will set all the mass scales equal to one. 

The action describing a (l + l)-dimensional gauge field interacting with point-like 
charges can be written as 

= J d 2 xF 2 + q J dafAp, (26) 

where F is the field strength and x M is the coordinate of a point charge q. With 
q =constant, the above action is gauge invariant under — > + d^u. However, 
we would like to promote the charge q to a function of a scalar field 0, that is, q = <f> . 
With this substitution, however, the action (|26|) is no longer gauge invariant. In 
order to restore the gauge symmetry, we shall modify it in the following way: 

S 1+1 = J d 2 xF 2 + J <f) N dafTl^A", (27) 

where U^ u = i]^ u — is the transverse projector. The modification is only 
significant for a varying 0. For a constant q = <p N , the actions (}2*Tj) and (j2Sj) are 
equivalent. 

The generalization of this action to a (3 + l)-dimensional model with 2-branes 
is straightforward: 

S 3+1 = J d 4 xF 2 + J <f> N da^U^'U^'U^' A^y (28) 

where da^ is the world- volume element. This action is manifestly gauge invariant 
under (jSJ). 

An unusual feature of the actions (|27|). (j2Hj) is that the projection operators 
are non-local. It appears that such operators cannot be obtained from a local field 
theory by integrating out a finite number of heavy particles. It Ref . j3| , the action 
(J27)l was obtained after integrating over massless fermions in a (1 + l)-dimensional 
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model. However, it is not clear how generic that model is and whether or not it can 
be extended to higher dimensions. These issues require further study. 

As we already mentioned, the branes (or charges) in question can be regarded 
as fundamental objects, or as solitons of the effective field theory. We shall briefly 
review the latter possibility. Following j3J, we assume that in the 4D low-energy 
effective theory F is mixed with a certain phase field a, 

i-oF. (29) 

The interaction (|29|) is invariant under the shift symmetry 

a -> a + 2vr, (30) 

as well as under gauge transformation of the three-form field A^ va . The equation of 
motion of the three-form field then demands that the variations of the field strength 
and a with respect to any particular coordinate must satisfy 

AF = — Aa (31) 
2n 

Thus, any stable solitonic configuration across which a changes by a finite amount, 
inevitably acts as a source for the three-form field. Such configurations do indeed 
exist. Since a is a phase, the potential of a must respect the shift symmetry ([30jh 
As long as this is the case, irrespective of the precise form of this potential, due to 
topological reasons there must be domain wall solutions across which a changes by 
Aa = 2tt. Then, according to (J3T)) . such walls acquire a charge q under A. Thus, 
the values of the field strength on the two sides of the wall differ by 

AF = q. (32) 

The charge q defines the minimal step by which F can change from one region of 
the Universe to another. 

For our purposes, however, q cannot be truly constant, because it has to depend 
on 0. Also, if we want to think of a as a phase of a certain complex order parameter, 

X = \X\e ia , (33) 

with a VEV around the Planck scale, |X| ~ 1, then q will depend on \X\ as well. 
In terms of the fields X and 0, the gauge-invariant coupling (J28j) can be written as 

A^ n""' IT"' IF 7 ' (ie^ya (Xd a X^ - X^d a X) <t> N ) (34) 

This is the lowest possible a — A- mixing operator invariant under gauge (jHJ) and shift 
()3U|) symmetries, and under the Z 2 n symmetry which acts on the field as 

-»• e*# . (35) 
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The action of this symmetry on the other fields can be defined in two alternative 
ways. One possibility is to demand 

a -»• -a, A -> A, (36) 

or equivalently, X — > X>. Note that the transformation a — > —a replaces solitons 
by anti-solitons and vice versa. An alternative choice would be to demand 

a -> a, A — > -A (37) 

Eq (|34j) is invariant under both of these choices. At low energies, where and 
|X| can be treated as constants, the above coupling reduces to the one of ()29|). 
Note that we are not demanding any approximate or exact U(l) symmetry under 
the shift a —>■ a+constant. For the existence of the wall, all that we need is that 
\X\ = be a maximum of the potential. This suffices to ensure the existence of a 
stable configuration across which a changes by 2ir. The corresponding change of F 
through the wall will be 

AF ~ (f) N . (38) 



B F-waves from the walls 

In 3 + 1-dimensions, gauge-invariant free 3-form fields (just as the (l + l)-dimensional 
electromagnetic field) have no propagating degrees of freedom, so there are no wave 
solution. The situation is different in our framework, where brane charges depend 
on the field 0. At the end of inflation, rolls away from zero, and the brane charges 
change in time. This change triggers the corresponding change of the three-form 
field strength, which propagates away from the brane in the form of a shock wave. 
We shall now discuss this dynamics. 

To illustrate the point, we shall restrict ourselves to a (1 + l)-dimensional ex- 
ample. Generalization to (3 + 1) dimensions is straightforward. Labelling the two 
spacetime coordinates x 11 by t and z, we shall consider an isolated static "brane" 
located at z = 0. The gauge- invariant Lagrangian of interest is 

L = - 1 -F 2 + A,W\<p N 5{z))+\d^ - (-ml ff \<f>\ 2 + |0| 4 ), (39) 

where m e ff is the effective mass of after inflation, and we have ignored higher- 
order terms in the potential of 0. When rolls away from = 0, the change of the 
brane charge induces a back reaction on <ft. This back reaction is suppressed by the 
brane charge ~ <f> . Since in the domains of interest, the final VEV of is very 
small, the back reaction on is negligible, and we shall ignore it. 

Thus, we shall study the dynamics of a gauge field in the background of a time- 
dependent charge. This dynamics is governed by the following equation 

0" F, v = rU (<f> N (t) S(z)) = Vuo 4> N (t) 5{z) - J e"*« <f N (po), (40) 
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where p is the two-momentum, and <Pn{Po) is the Fourier-transform of <ft N (t). A 
straightforward integration gives the following two equations 

d z F w = l -d z (<p N (t - z)6(z) -<f(t + z)6(-z)), (41) 

d t F w = ld t (<t> N {t-z)9{z)-<t> N {t + z)9{-z)), (42) 
which are solved by 

F w = \ ( <j> N (t - z) 6(z) - <P N (t + z) 9{-z) ) . (43) 

This solution describes waves propagating away from the brane in two opposite 
directions. From this solution it is also obvious that after the transition is finished 
and assumes a constant value (fitoday, the change of F across the brane is given by 

AF = <f>Zday (44) 
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